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Abstract 

We calculate the cosmic ray positron and antiproton spectra of a gauge-Higgs dark matter 
candidate in a warped five-dimensional 50(5) x C/(l) gauge-Higgs unification model. The stability 



Qh' of the gauge-Higgs boson is guaranteed by the H parity under which only the Higgs boson is odd 

T^, ' at low energy. The 4-point vertices of HHW~^W~ and HHZZ, allowed by H parity conservation, 

have the same magnitude as in the standard model, which yields efficient annihilation rate for 

>. ... ^ _ 

CN niH > rnw- The most dominant annihilation channel is HH — )• W^W followed by the subsequent 

00 ' 

g : decays of the W bosons into positrons or quarks, which undergo fragmentation into antiproton. 

t^^ ■ Comparing with the observed positron and antiproton spectra with the PAMALA and Fermi/LAT, 

o. 

^^ , we found that the Higgs boson mass cannot be larger than 90 GeV, in order not to overrun the 

observations. Together with the constraint on not overdosing the Universe, the valid range of the 
^ ■ dark matter mass is restricted to 70 — 90 GeV. 



I. INTRODUCTION 



Since Lee and Weinberg put bound on a stable neutrino mass from the relic density of the 
Universe in 1977 [l|, various cosmological and astrophysical observations have inspired and 
also constrained many theoretical models in particle physics. In particular the very precise 
measurement of the cosmic microwave background radiation in the Wilkinson Microwave 
Anisotropy Probe (WMAP) experiment established the presence of cold dark matter (CDM) 
in our Universe |2|] . These data are urging new physics (NP) beyond the standard model 
(SM) to provide the stable CDM particle over cosmological time scale. 

Recently, a number of high energy cosmic ray experiments suggested the possibility of 
indirect detections of CDM annihilation in the galactic halo. The Payload for Antimatter 
Matter Exploration and Light-nuclei Astrophysics (PAMELA) collaboration reported the 
excess of positrons in its energy spectrum of 10 — 100 GeV over the expected high energy 
cosmic ray interaction with the interstellar medium [3|] . More mysterious result was that 
the antiproton data, on the other hand, seem to be consistent with the expected spectrum 
from the collisions between high-energy protons in cosmic rays and the nuclei of hydrogen 
and helium atoms in the interstellar medium j^ . In addition, the Fermi/LAT [5| and the 
HESS p] collaborations also reported a smooth but harder electron-flux spectrum than the 
expected background at 100-1000 GeV. 

In the literature there have been great efforts to explain these anomalous data by the 
annihilation of CDM particles in NP models, such as the lightest supersymmetric particle [7|], 
the Kaluza-Klein particle |8|, and other CDM candidates J9|]. It was shown that many CDM 
particle candidates such as in supersymmetric (SUSY) models and the universal extra di- 
mensional (UED) model require a larg e boost factor to explain the observed signal. Even 



u 



though the inhomogeneity of CDM |10| or the Sommerfeld mechanism [ill] can help, the val 



ues of these boost factors are limited. Some other NP models^!, e.g., a new long range 

n 

force in the dark sector pj| , can explain the PAMELA data without extreme boost factor. 
In this paper we adopt a different approach to these astrophysical data, i.e., using them to 
constrain a NP theory, instead of explaining them to support a theory. The positron excess 
could be explained by astrophysical sources like pulsars 13| or supernova remnants [ij]. 
Historically the astrophysical observations have played the role of a strong constraint on a 
NP model. For instance, the very precise WMAP measurement of the relic density eliminates 



a large portion of the parameter space of the constrained minimal supersymmetric SM to 



avoid overdosing the Universe 15] 



We will study the positron and antiproton energy spectra from CDM annihilation in 
a recently proposed 5*0(5) x f/(l) gauge- Higgs unification model based on a warped five- 



dimensional (5D) spacetime 16l-ll8|. Here the dark matter is nothing but the Higgs bo- 



son 



19| . which is a part of the fifth component of a gauge boson field in the model. The H 
parity, under which only the Higgs boson has odd parity at low energy, is preserved from 
the gauge structure of the theory. Triple vertices such as WWH, ZZH, and ffH (here / 



is a SM fermion) vanish to all orders in perturbation theory 19|, |20|, thus the Higgs boson is 
stable and becomes a dark matter candidate. Instead, the 4-point vertices of HHW~^W~, 
HHZZ, and HHff are allowed by the H parity conservation. The Higgs boson can be 
thermally produced in the early Universe via WW, ZZ, ff — ?■ HH, of which the rate is 
highly predictive with essentially one free parameter, the Higgs boson mass rriH- 

In Ref. 191], it was shown that the Higgs boson mass needs to be at 70 GeV in order to 
explain the WMAP data. When rriH < 70 GeV, the total annihilation rate is very small 
because the kinematically allowed 2 — )> 2 processes into light fermions are suppressed by 
small Yukawa couplings. The corresponding relic density becomes too big, which overdoses 
the Universe. If rriH > 70 GeV, the cross section of HH -^ WW^*^ becomes very large, 
since the magnitude of HHWW coupling is the same as in the SM. The corresponding relic 
density is too low. Nevertheless, this is not ruled out. It is possible that the dark matter 
can be produced nonthermally from, e.g., other quasi-stable Kaluza-Klein states or other 
topological objects. In this work, we assume that the thermal source is not the only source 
of relic dark matter, and so Higgs boson mass can be larger than 70 GeV. The present upper 
bound on the Higgs mass comes from the unitarity requirement, which limits rriH < 0{1) 



TeV 



m. 



We note that this large Higgs mass yields too efficient annihilation into W^W~ and 
possibly ZZ, which leads to potentially large positron and antiproton signals. Using the 
observed positron and antiproton spectra in the PAMELA and Fermi/LAT experiments, we 
can set the upper limit on the Higgs boson mass. On the other hand, the collider signal for 



this model is way too small for detection 22]. 



In this work, we study the positron and antiproton spectra from annihilation of the 



gauge-Higgs dark matter in the halo. We use the cosmic ray propagation code Galprop 23] 



to calculate the propagation of the positron and antiproton from the halo to the Earth, and 
compare with the spectra measured by PAMELA 3|, \^. The most dominant annihilation 
channel for niH > 70 GeV is HH — ?> W~^W~, followed by the subsequent decays of the W 
bosons into positrons or quarks. The quarks undergo fragmentation into antiproton. Here 
we do not attempt to explain the anomaly observed by PAMELA. Instead we use them as 
the constraints on the model: the resulting spectra obtained from the CDM annihilation 
should not exceed the ones measured by PAMELA. 

It is well known that the process HH — t- W^W~ grows with the center-of-mass (cm.) 
energy -\/i. Naively, the longitudinal polarization of the W boson behaves like p^/rriw when 
y/s ^ mw Therefore, we expect the scattering amplitude squared grows as s'^/m^r before 
any UV physics comes in to unitarize the theory.^ For the CDM annihilation in the halo 
where v ~ 10"'^, the cm. energy is just roughly 2mH- Thus, we expect the annihilation 
rate will grow with the Higgs boson mass, as rapidly as rn^j/rn^-, so long as mu < 1 TeV. 
We will show that when the Higgs boson mass is 90 GeV or above, the resulting positron 
spectrum is already well above that measured by PAMELA. Similar conclusion holds true 
for the antiproton spectrum. We can, therefore, conclude that the Higgs boson mass cannot 
be larger than about 90 GeV in this gauge-Higgs model. This is the main result of our work. 

The organization of the paper is as follows. In the next section, we briefly review the 
model and describe the effective interactions used in this calculation. We give details about 
the calculation of the positron spectrum and antiproton spectrum in Secs. lIIII and IIVI re- 
spectively. We also perform the comparison with the measured spectra. We conclude in 
SeclVl 

II. THE EFFECTIVE INTERACTIONS IN THE 50(5) x [/(I) GAUGE-HIGGS 
MODEL 

The model under consideration is a 5*0(5) x f/(l) gauge-Higgs unification model in the 
5D Randall- Sundrum warped space [l6|. The Higgs boson is the fluctuation mode of the 
Aharonov-Bohm phase 9h along the fifth dimension 



2^, i.e., Oh = eH + H{x)/fH. The 



^ In this gaugc-Higgs model based on five-dimensional spacetime, it is the Kaluza-Klein states (of order 
TeV) of the gauge bosons which unitarize the scattering amplitude. 



four- dimensional effective Lagrangian of the Higgs boson is 

C = V,s{Oh) - ml,{eH)W;W~^ - \nl{eH)Z^,Z^^ -J2mfi9H)i^f^f , (1) 

/ 

where the mass functions are 

mwiOn) = -gfnsineH, mziOn) = -gzfHsm9H, rnf{9H) = yfiusmdH ■ (2) 

Here g is the weak gauge couphng and gz = g/ cos 9w 

The effective potential Vcs{6h) of the Higgs boson is generated at one loop level. It is 



finite and cutoff independent. As shown in Ref. 16|, the large contribution of 5D top quark 
field changes the global minimum of VesiOn) into 9h = ±vr/2: the W, Z gauge bosons as well 
as the SM fermions acquire their masses (see Eq.([2])) and thus the electroweak symmetry is 
broken dynamically. 

In this model the global minimum at 6h = ±7r/2 dynamically generates a new H parity, 
under which the Higgs boson has odd parity while all the other SM particles have even 
parity. This iZ-parity prohibits triple vertices of the Higgs boson with the SM particles, 
such as HW~^W~ , HZZ and Hff: it preserves the stability of the Higgs boson so that the 
Higgs boson can be a CDM candidate. 

At low energy this model has two parameters, fn and rriH- The parameter fn is deter- 
mined by the observed mw and mz, i.e., fn ~ 246 GeV. The value of mn is, in principle, 
determined if the whole matter content in the model is fixed in detail. Without a priori 
knowledge of UV physics, we treat mn as a free parameter. 

Because of the absence of the triple vertices of the Higgs boson with the SM particles, we 
do not have significant constraint on rriH from collider physics phenomenology. Instead the 
observed relic density of CDM in the WMAP experiment can provide a meaningful one [19|] . 
If the Higgs boson mass becomes heavier, their annihilation into W^W~ and possibly into 
ZZ are kinematically accessible, leading to smaller relic density of the Universe. On the 
other hand, lighter mu yields too small annihilation cross section, which is excluded as 
overdosing the Universe. The Higgs boson mass mu = 70 GeV can explain the observed 
relic density. We employ the WMAP data as not overdosing the universe: the Higgs boson 
mass can be heavier than 70 GeV. 

The effective interactions used in this work are 

£ = \glH^Z,Z^ + \gVw;W-^ + J^ |^i/^^/V^/. (3) 



Note that the HHW^W~ and HHZZ vertices have the same couphngs with the SM ones 
except for the opposite sign. Since the annihilation cross section of HH — )■ W~^W~ , ZZ 
increases with the cm. energy -y/i (or mn when the relative velocity of two if 's is very 
small), heavier Higgs boson could leave too much excess of positrons and antiprotons from 
W or Z decays. 

III. POSITRON SPECTRUM 

The dominant process for the cosmic ray positrons from the Higgs boson annihilation is 
the leptonic decay of the W~^ boson: 

HH -^ W+W- ->e+ + u, + X. (4) 

The next dominant process is HH — > ZZ — )■ e~^ + X . Positrons can also come from 
the hadrons, which are the fragmentation products of the quarks from W decays. These 
positrons are much softer than those coming directly from the W boson decay. We shall ig- 
nore these soft positrons. There are also other processes HH — )■ // (for f = b,c) to produce 
positrons in the subsequent decays of the fragmentation products, but they are certainly 
subleading because of the small Yukawa couplings compared with the gauge coupling. 
A back-of-envelope calculation shows that the annihilation rate is given by 

<-">— - "^(^^ - ^^) (2^") = it (=» - 4 ^ 4) ■ <'' 

where 2/3// = 2a/1 — Am\j/ s is the relative velocity of the two non-relativistic incoming 
Higgs bosons in their cm. frame, and I3w,z = \ h- ~ ^f^w zl ^- ^^ i^ easy to see that the 
annihilation rate grows as s/m"^^. 

The positron flux observed at the Earth is given by 

^AE) = ^-^feAE), (7) 

with Ve+ is close to the velocity of light c. The function /e+ (E) satisfies the diffusion equation 
of 

^ - KiE)V'f - ^ ibiE)f) = g , (8) 
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where the diffusion coefficient is K{E) = Ko{E/GeVY and the energy loss coefficient is 
b{E) = E'^/{GeV x te) with te = 10^^ sec. The source term Q due to the annihilation is 

where rj = 1/2 for identical scalar DM particle in the initial state. The summation is over 
all possible channels that can produce positrons in the final state, and dNf,+ /dEf,+ denotes 
the spectrum of the positron energy per annihilation in that particular channel. 
In our analysis, the source term is given by 



^ 1 / PCDM 



2 VM, 



CDM 



{av)HH^WW —TT=, y- {<yv)HH~,ZZ 



dE.+ ' /....-.^^ ^^ 



(10) 



where the normalizations of N^^ and N^+^ are 



e^ 



/ ^—^dx = B{W+ -^ e+z/e), / ^-f^dx = 2xB{Z^ e+e"). (11) 

We first calculate the energy spectrum of the positron in the W~^ or Z rest frame, then 
boost it to the cm. frame of the HH system. In the calculation we include all the off-shell 
effects of the W~^ and W~ bosons and in both Z bosons. We then put the source term into 

n 

Galprop [23j] to solve the diffusion equation. 

Note that the contribution from HH — )■ ZZ channel is subleading in our calculation. It 
only accounts for about 1% and 4% of the contribution from HH — )■ W'^W~ channel for 
rriH = 82 and 90 GeV, respectively. For mn = 100 GeV, the ZZ channel contribution is 
as large as 20% of the WW one. It is easy to understand: for rriH below 90 GeV the ZZ 
channel is below the production threshold and thus one of the Z bosons has to be off-shell, 
while at rriH = 100 GeV both Z bosons are already on shell. 

Figure [1] shows the positron spectra for rriH = 70, 82, 90, and 100 GeV with the measured 
spectrum of PAMELA [31]. For rriH = 70 GeV, the gauge-Higgs CDM annihilation cannot 
explain the rising feature in the positron spectrum. The annihilation rate is too low since the 
W^W~ channel is not kinematically open yet. When the Higgs boson mass is above mw, 
the annihilation rate increases quickly and can explain part of the rise-up in the spectrum, 
as shown by the rriH = 82 GeV curve. Yet, it still cannot explain the two highest energy 
points because the mass of the Higgs boson is not heavy enough. Once the Higgs boson 
mass rises to 90 — 100 GeV, it passes through the second highest energy point, but the 
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FIG. 1: The positron- fraction spectrum predicted for the gauge-Higgs dark matter model for 
rriH = 70, 82, 90, 100 GeV, assuming that it accounts for all the dark matter of the Universe. The 
PAMELA data are also shown. 



annihilation rate increases far more than the lower part of the spectrum allowed. Just by 
visual checking, rriH larger than 90 GeV is easily ruled out by the measured spectrum. 

We also show the total electron and positron flux in Fig. [21 Our results are compared 
with the measurement of the cosmic ray e"*" + e~ spectra from 20 GeV to 1 TeV with the 
Fermi LAT |5|, and from 340 GeV with the HESS p]. These cosmic ray electron spectra are 
distributed over very high energy up to 1 TeV. Unless the CDM mass is very heavy about 
800 — 1000 GeV, the whole energy spectrum cannot be explained. Instead of explaining 
high energy spectrum, we focus on the region of 20 GeV < E < 100 GeV. The case of 
irtii = 82 GeV is marginally allowed by the data while the rriH = 90 GeV case already 
outruns the observation. 
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FIG. 2: Total flux of electron and positron predicted for the gauge-Higgs dark matter model for 
mn = 70, 82, 90, 100 GeV. Data from Fermi/LAT and HESS are shown. 

IV. ANTIPROTON SPECTRUM 



Similar to the treatment for positron flux, the antiproton flux can be obtained by solving 
the diffusion equation with an appropriate source term for the input antiproton spectrum: 



Qa 



V 



Pdn 



>r^ dNp 



dTfi 



(12) 



where rj = 1/2, and Tp is the kinetic energy of the antiproton which is conventionally used 
instead of the total energy. We solve the diffusion equation using Galprop [23 1. 

In our case, the dominant contribution to the cosmic ray antiproton production comes 
from 

HH ^W+W- ^ {qq') {qq') -^ p + X , (13) 

followed by HH — )■ ZZ — )■ {qq){q'q') — )■ p + X. In the last step, we adopt a publicly available 
code 251] to calculate the fragmentation function Dg_^h{z) for any quark q into hadrons h, 
e.g., p,p,7r. The fragmentation function is then convoluted with energy spectrum dN/dE 
of the light quark to obtain the energy spectrum of the antiproton dN/dEp. The next 
contribution comes from HH -^ bb ^ p + X. Since the annihilation is smaller by two 
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FIG. 3: The fraction of the antiproton p/{p + p) versus the kinetic energy of the antiproton 
predicted for the gauge-Higgs dark matter model for rriH = 82, 90, 100 GeV. The antiproton data 
of PAMELA are also shown. 

orders of magnitude, we ignore this and the other subleading contributions. The source 
term dN/dTp is then implemented into Galprop to calculate the propagation from the halo 
to the Earth. 

The resulting antiproton fraction as a function of its energy is shown in Fig. [3l Three 
theoretical curves are for niH = 82, 90, 100 GeV, compared with the observed PAMELA 
data. Note that in this figure we only include the leading contribution from HH — > W^W~ 
because the next subleading contribution only accounts for less than 5% of the leading one 
for rriH < 90 GeV. We can see that the rriH = 82 (just above the WW threshold) curve is 
barely consistent with the data. On the other hand, the rriH = 90 and 100 GeV curves are 
obviously over the data. 

Since the theoretical background curve contains relatively high uncertainties due to dif- 
ferent models used in Galprop, we do not perform any confidence- level exclusion analysis. 
Instead, we only use a visual inspection. Very conservatively, we conclude that the rriH = 90 
GeV or above is excluded by the antiproton spectrum. A similar conclusion can be drawn 
from the positron spectrum. Overall, rriH = 90 GeV or above is strongly disfavored by the 
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data 



V. CONCLUSIONS 

In the 5*0(5) x f/(l) gauge-Higgs unification model based on the Randall-Sundrum warped 
spacetime, the Higgs boson is a cold dark matter candidate. Its stability over the cosmolog- 
ical time scale is guaranteed by the H parity under which only the Higgs boson has negative 
parity at low energy. The triple vertices of HWW, HZZ and Hff vanish. However the 
4-point vertices of HHWW, HHZZ, and HHff are allowed. 

Interpreting the WMAP data on the relic density as not overdosing the Universe, the 
Higgs boson mass is constrained from below as niH > 70 GeV. Noting that the heavy Higgs 
boson generates more efficient annihilation into W^W~ and ZZ, we study the cosmic ray 
positron and antiproton spectra and compare them with the PAMELA and Fermi/LAT 
observations. We do not aim at explaining the up-rising positron spectrum observed with 
the PAMELA. Instead we use the data to constrain the model. The Higgs boson annihilation 
rate is shown to increase with increasing Higgs boson mass. The mn = 82 GeV is marginally 
consistent with both the observed positron and antiproton spectra. However, the Higgs 
boson mass of 90 GeV or more is obviously ruled out by the observed data. 

Our purpose is to constrain the 5*0(5) x f/(l) gauge-Higgs unification model by the 
PAMELA and Fermi/LAT data, which turned out to be very significant. The present limit 
on the model comes from the unitarity requirement, which limits the Higgs boson to be 



less than about 0{1) TeV 21|. We show explicitly in this paper that the Higgs boson mass 
cannot be larger than 90 GeV, otherwise it overruns the observed positron and antiproton 
spectra of PAMELA and Fermi/LAT. Together with the constraint on the closure of the 
Universe, the gauge-Higgs boson mass as CDM is now restricted to 70 — 90 GeV. 
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